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rality of values each representing a condition of the re- 
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to correct extraneous values in the data which could ad- 
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Description 

Cross-Reference To Related Patent and Application 

5 [0001] Related subject matter is disclosed in a copending U.S. patent application of Harry Yang, Daniel L. Schwarz, 
Christopher M. Embres, Richard L. Moore. Penry D. Haaland and Paula V. Johnson entitled "Computerized Method 
and Apparatus for Analyzing Readings of Nucleic Acid Assays". Serial No. 09/196.123. filed on November 20, 1998, 
and in U.S. Patent No. 6.043.880 of Jeffrey P. Andrews, Christian V. 0*Keefe, Brian G. Scrivens, Willard C. Pope, 
Timothy Hansen and Frank L. Failing entitled "Automated Optical Reader for Nucleic Acid Assays", the entire contents 

10 of said application and patent being expressly incorporated herein by reference. 

BACKGROUND OF THE INVENTION 
Field of the Invention : 

15 

[0002] The present invention relates generally to a computerized method and apparatus for analyzing sets of readings 
taken of respective samples in a biological or chemical assay, such as a nucleic acid assay, to determine which samples 
possess a certain predetermined characteristic. More particularly, the present Invention relates to a computerized meth- 
od and apparatus which graphically plots optical readings of a biological or chemical sample taken at different times 
20 during a reading period, corrects for an additive background value present in the readings, and compares the corrected 
readings to a threshold to detect for the presence of a particular characteristic of the sample, such as the presence of 
a targeted pathogen in the sample. 

Description of the Related Art : 

25 

[0003] In the clinical diagnosis of bacterial diseases, including sexually transmitted diseases such as gonorrhea (GC) 
and chlamydia trachomatis (CT), a sample of body fluid obtained from the patient is cultured to test for the presence 
of the particular bacterium of interest. Unfortunately, this is a relatively time-consuming process, generally requiring 
several days to produce a definitive result. During this time, a patient suspected of having such a disease must be 

30 isolated to prevent further spread of the disease. 

[0004] The advent of DNA probes, which can identify specific bacteria by testing for the presence of a unique bacterial 
DNA sequence in the sample obtained from the patient, has greatly increased the speed and reliability of clinical 
diagnostic testing. A test for the presence of CT or GC, for example, in a sample can be completed within an hour or 
less using DNA probe technology. This allows treatment to begin more quickly, and avoids the need for long patient 

35 isolation times. 

[0005] In the use of DNA probes for clinical diagnostic purposes, a nucleic acid amplification reaction is usually 
carried out in order to multiply the target nucleic acid into many copies or ampltcons. Examples of nucleic acid ampli- 
fication reactions Include strand displacement amplification (SDA) and polymerase chain reaction (PCR). Unlike PCR, 
SDA is an isothermal process which does not require any external control over the progress of the reaction which 
40 causes amplification. Detection of the nucleic acid amplicons can be carried out in several ways, all involving hybrid- 
ization (binding) between the target DNA and specific probes. 

Many common DNA probe detection methods involve the use of fluorescein dyes. One known detection method is 
fluorescein energy transfer. In this method, a detector probe is labeled both with a fluorescein dye that emits light when 
excited by an outside source, and with a quencher which suppresses the emission of light from the flourescein dye in 
45 its native state. When DNA amplicons are present, the fluourescein-labled probe binds to the amplicons, is extended, 
and allows fluorescent emission to occur. The increase of fluorescence is taken as an indication that the disease- 
causing bacterium is present in the patient sample. 

[0006] Several types of optical readers or scanners exist which are capable of exciting fluids samples with light, and 
then detecting any light that is generated by the fluid samples in response to the excitation. For example, an X-Y plate 

50 scanning apparatus, such as the CytoFluor Series 4000 made by PerSeptive Biosystems, is capable of scanning a 
plurality of fluid samples stored in an array of microwells. The apparatus includes a scanning head for emitting light 
towards a particular sample, and for detecting light generated from that sample. During operation, the optical head is 
moved to a suitable position with respect to one of the sample wells. A light emitting device is activated to transmit 
light through the optical head toward the sample well. If the fluid sample In the well fluoresces in response to the emitted 

55 light, the fluorescent light is received by the scanning head and transmitted to an optical detector. The detected light 
is converted by the optical detector into an electrical signal, the magnitude of which is indicative of the intensity of the 
detected light. This electrical signal is processed by a computer to determine whether the target DNA is present or 
absent in the fluid sample based on the magnitude of the electrical signal. Each well in the microwell tray (e.g.. 96 



2 




EP1 158 449 A2 

microwells total) can be read in this manner. 

[0007] Anothier more efficient and versatile sample well reading apparatus known as the BDProbeTec® ET system 
manufactured by Becton Dicl<inson and Company is described in the above-referenced U.S. Patent No. 6,043,880. In 
that system, a mlcrowell array, such as the standard microwell array having 12 columns of eight microwells each (96 
5 microwells total), is placed in a moveable stage which is driven past a scanning bar. The scanning bar includes eight 
light emitting/detecting ports that are spaced from each other at a distance substantially corresponding to the distance 
at which the microwells in each column are spaced from each other. Hence, an entire column of sample microwells 
can be read with each movement of the stage. 

[0008] As described in more detail below, the stage is moved back and forth over the light sensing bar, so that a 
10 plurality of readings of each sample microwell are taken at desired intervals. In one example, readings of each microwell 
are taken at one minute intervals for a period of one hour. Accordingly, 60 readings of each microwell are taken during 
a well reading period. These readings are then used to determine which samples contain the particular targeted disease 

or diseases (e.g., CT and/or GC). 

[0009] Several methods are known for analyzing the sample well reading data to determine whether a sample con* 
15 tained in the sample well includes the targeted disease or diseases. For instance, as discussed above, a nucleic acid 
amplification reaction will cause the target nucleic acid (e.g., CTorGC) to multiply into many amplicons. Thefluorescein- 
labeled probe which binds to the amplicons will fluoresce when excited with light. As the number of amplicons increases 
overtime while the nucleic acid amplification reaction progresses, the amount of fluorescence correspondingly increas- 
es. Accordingly, after a predetermined period of time has elapsed (e.g. 1 hour), the magnitude of fluorescence emission 
20 from a sample having the targeted disease (a positive sample) is much greater then the magnitude of fluorescence 
emission from a sample not having the targeted disease (a negative sample). In actuality, the magnitude of fluorescence 
of a negative sample essentially does not change throughout the duration of the test. 

[001 0] Therefore, the value of the last reading taken for each sample can be compared with a known threshold value, 
which has previously been determined. If the sample value is above the threshold value, the sample is identified as a 
25 positive sample in which the targeted disease is present. However, if the last value taken of the sample is below the 
threshold value, the sample is identified as a negative sample free from the disease. 

[001 1] Although this "endpoint detection" method can generally be effective in identifying positive and negative sam- 
ples, it is not uncommon for this method to incorrectly identify a negative sample as being positive or vice versa. That 
is, the accuracy of the value of any individual sample reading can be adversely effected by factors such as a bubble 

30 forming in the sample, obstruction of excitation light and/or fluorescence emission firom the sample due to the presence 
of debris on the optical reader, and so on. Accordingly, if the final reading of a particular sample is erroneous and only 
that reading is analyzed, the likelihood of obtaining a false positive or false negative result is high. 
[0012] Furthermore, in some instances, the amount of amplicons decreases as the testing progresses. Due to this 
decrease in amplicons, the magnitude of fluorescent emission from the sample decreases with the passage of time. 

35 Accordingly, the magnitude of the last reading taken of the sample can be less than the magnitude of a reading taken 
at the time when the amount of amplicons present in the sample is at its peak. In some instances, the decrease in 
amplicons can result in the magnitude of fluorescent emission from the sample being lower then the predetermined 
threshold value, in which event the positive sample is falsely identified as a negative sample. 
[0013] In order to avoid these drawbacks, other methods have been developed. In one method, the overall change 

40 in the magnitudes of sample readings is calculated and compared to a known value having a magnitude indicative of 
a positive result. Accordingly, if the magnitude of change is greater than the predetermined value, the sample is iden- 
tified as a positive sample having the targeted disease. On the other hand, if the magnitude of change is less than the 
predetermined value, the sample is identified as a negative sample. 

[0014] Although this method may be more effective then the endpoint detection method discussed above, certain 
45 flaws in this method also exist. For example, if a sample contains a particularly large amount of the disease, the amount 

of amplicons generated due to the amplification process may reach a maximum at the time the initial reading is taken. 

and increase very little, if at all, or decrease, throughout the duration of the reading period. In this event, the change 

which occurs between the initial readings and final readings is minimal even though the sample is very positive. Hence, 

the sample may incorrectly be identified as a negative sample, 
50 [001 5] Another known method is the acceleration index method which measures incremental changes in the sample 

readings, and compares those changes to a predetermined value. Although this method is generally effective, the 

accuracy of its results are susceptible to errors present in the individual readings, as well as readings taken of samples 

in which the amount of amplicons decreases over time. 

[0016] Accordingly, a continuing need exist for a method an apparatus for analyzing data representative of readings 
55 taken of sample wells to accurately identify the samples as being positive or negative for a particular disease. 
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SUMMARY OF THE INVENTION 

[0017] An object of the present invention is to provide a method and apparatus for accurately interpreting the values 
of data obtained from taking readings of a biological or chemical sample, to ascertain the presence of a particular 

5 disease in the sample based on the data values. 

[0018] Another object of the invention is to provide a method and apparatus, for use with an optical sample well 
reader, which accurately interprets data representing magnitudes of fluorescence emissions detected from the sample 
at predetermined periods of time, to ascertain the presence of a particular disease in the sample. 
[001 9] A further object of the invention is to provide a method and apparatus for analyzing data obtained from reading 

10 a biological or chemical sample contained in a sample well, and without using complicated arithmetic computations, 
correcting for errors in the data which could adversely affect the results of the analysis. 

[0020] These and other objects of the invention are substantially achieved by providing a computerized method and 
apparatus for analyzing numerical data pertaining to a sample assay comprising at least one biological or chemical 
sample, with the data including a set of data pertaining to each respective sample, and each set of data including a 

15 plurality of values each representing a condition of the respective sample at a point in time. The method and apparatus 
assigns a respective numerical value to each of the data values, removes an additive background value from each of 
the data values to produce corrected data values, compares the values to a threshold value to determine when the 
data values begin to exceed the threshold value, and controls the system to indicate whether the sample has a pre- 
determined characteristic based on a result of the comparison. Additionally, prior to the comparison to the threshold 

20 value, filtering, normalizing and other correcting operations can be performed on the data to correct extraneous values 
in the data which could adversely affect the accuracy of the results. 

[0021 ] The method and apparatus can perform all of the above functions by representing each of the plurality of data 
values as points on a graph having a vertical axis representing the magnitudes of the values and a horizontal axis 
representing a period of time during which readings of the sample were taken to obtain said plurality of data values, 
25 correcting the data values to eliminate an additive background value present in each of the data values to produce a 
corrected plot of points on the graph, with each of the points of the corrected plot of points representing a magnitude 
of a corresponding one of the values. The corrected plot of points is compared to a threshold value to determine if and 
when the points exceed the threshold value, to therefore determine whether a certain condition, such as a pathogen, 
exists in the sample to which the set of data pertains. 

30 

BRIEF DESCRIPTION OF THE DRAWINGS 



[0022] The and other objects of the invention will be more readily appreciated from the following detailed description 
when read in conjunction with the accompanying drawings in which: 



35 



Fig. 1 is schematic view of an apparatus for optically reading sample wells of a sample well array, and which 
employs an embodiment of the present invention to interpret the sample well readings; 

Fig. 2 is an exploded perspective view of a sample well tray for use in the sample well reading apparatus shown 
in Fig. 1; 

40 Fig. 3 is a detailed perspective view of a stage assembly employed in the apparatus shown in Fig. 1 , for receiving 

and conveying a sample well tray assembly shown in Fig. 2. 

Fig. 4 is a diagram illustrating the layout of a light sensor bar and corresponding fiber optic cables, light emitting 
diodes and light detector employed in the apparatus shown in Fig. 1 , in relation to a sample well tray being conveyed 
past the light sensor bar by the stage assembly shown in Fig. 3; 
45 Fig. 5 is a graph illustrating values representing the magnitudes of fluorescent emissions detected from a sample 

well of the sample welt tray shown in Fig. 2 by the apparatus shown in Fig. 1 , with the values being plotted as a 
function of the times at which their corresponding fluorescent emissions were detected; 

Fig. 6 is a flowchart showing steps of a method for normalizing, filtering, adjusting and interpreting the data in the 

graph shown in Fig. 5 according to an embodiment of the present invention; 
50 Fig. 7 is a flowchart showing steps of the dark correction processing step of the flowchart shown in Fig. 6; 

Fig. 8 is a flowchart showing steps of the impulse noise filter processing step of the flowchart shown in Fig. 6; 

Fig. 9 is a flowchart showing steps of the dynamic normalization processing step of the flowchart show in Fig. 6; 

Fig. 10 is a flowchart showing steps of the impulse noise filter processing step of the flowchart shown in Fig. 6; 

Fig. 11 Is a graph which results after performing the dark correction, impulse noise filter, and dynamic normalization 
55 steps in the flowchart shown in Fig. 6 on the graph shown in Fig. 5; 

Fig. 1 2 is a flowchart showing steps of the step location and removal processing step of the flowchart shown in Fig. 6; 

Fig. 13 is a graph which results from performing the step location and repair step of the flowchart show in Fig. 6 

on the graph shown in Fig. 11; 
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Fig. 14 is a flowchart showing steps of the periodic noise filter operation step of the flowchart shown in Fig. 6; 
Fig. 15 is a flowchart showing steps of the well present detemnining step of the flowchart shown in Fig. 6; 
Fig. 16 Is a flowchart showing steps of the background correction step of the flowchart shown in Fig. 6; 
Fig. 17 is a graph which results from performing the background correction step of the flowchart show in Fig. 6 on 
the graph show in Fig. 13; 

Fig. 18 is a flowchart showing steps of the threshold comparison step of the flowchart shown in Fig. 6; 
Fig. 19 illustrates the comparison of a threshold value to the graph shown in Fig. 17; 

Fig. 20 illustrates the comparison of a threshold value to a graph generated from a very high positive sample; and 
Fig. 21 illustrates the comparison of a threshold value to a graph generated from a negative sample. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 



[0023] A well reading apparatus 100 according to an embodiment of the present invention is shown in Fig. 1. The 
apparatus 100 includes a key pad 102, which enables an operator to enter data and thus control operation of the 
apparatus 100. The apparatus 100 further includes a display screen 104, such as an LCD display screen or the like, 
for displaying "soft keys" which allow the operator to enter data and control operation of the apparatus 100, and for 
displaying information in response to the operator's commands, as well as data pertaining to the scanning information 
gathered from the samples in the manner described below. The apparatus also includes a disk drive 106 into which 
can be inserted a floppy disk for storing data generated by the apparatus 1 00, or from which the apparatus can read data. 

20 [0024] The apparatus 1 00 further includes as door 1 08 which allows access to a stage assembly 1 1 0 into which can 
be loaded a sample tray assembly 112. As shown in Fig. 2, a sample tray assembly 112 includes a tray 114 into which 
is loaded a microwell array 116, which can be a standard mtcrowell array having 96 individual microwells 118 arranged 
in 12 columns of 8 microwells each. The tray 114 has openings 120 which pass entirely through the tray and are 
arranged in 12 columns of eight microwells each, such that each opening 120 accommodates a microwell 118 of 

25 microwell array 116. After the samples have been placed into the microwells 118, a cover 122 can be secured over 
microwells 118 to retain each fluid sample in its respective microwell 118. Further details of the sample tray assembly 
112 and of sample collection techniques are described in the aforementioned U.S. Patent No. 6,043,880. 
[0025] Each microwell can include one type of detector probe, as described above, for identifying a particular disease 
(e.g., GC or CT). If the microwell array 11 6 is to be used to test for GC and CT in each patient sample, the microwells 

30 118 are arranged in groups of three microwells each, with one microwell in the group containing a reagent which is 
used to identify the presence of GC, another microwell in the group a reagent which is used to identify the presence 
of CT, and the third microwell containing an amplification control reagent AC, the purpose of which is described in more 
detail below. A fluid sample from a particular patient is placed in all three wells of a particular group of wells. 
[0026] Additionally, some (e.g., two) of the 96 microwells 118 in the microwell array 116 can be designated as control 

35 sample wells for a particular disease, such as CT, with one of the control sample wells containing a control positive 
sample and the other control well containing a control negative sample, the purpose of which is described In detail 
below. Also, other (e.g., two other) microwells 118 can be designated as control microwell sample wells for GC, with 
one of the control microwell sample wells including a positive sample and the other including a negative sample. Fur- 
thermore, other (e.g., two other) microwells 118 can be designated as control microwell sample wells for AC, with one 

40 of the control microwell sample wells including a positive sample and the other including a negative sample. Accordingly, 
in this example, a maximum of 30 patient samples can be tested for each microwell array 116 arranged in this manner 
(i.e., 30 samples x 3 microwells each = 90 microwells, with the remaining six microwells being used for control samples 
as discussed above). 

[0027] After the patient fluid samples have been placed in the appropriate microwells 118 of the microwell array 116 
45 in sample tray assembly 112, the sample tray assembly 112 is loaded into the stage assembly 110 of the well reading 
apparatus 1 00. The stage assembly 110 is shown in more detail in Fig. 3. Specifically the stage assembly 110 includes 
an opening 1 24 for receiving a sample tray assembly 1 1 2. The stage assembly 110 further Includes a plurality of control 
wells 126 which are used in calibrating and verifying the integrity of the reading components of the well reading appa- 
ratus 100. Among these control wells 126 is a column of eight normalization wells 127, the purpose of which is described 
50 in more detail below. The stage assembly 110 further includes a cover 128 which covers the sample tray assembly 
112 and control wells 126 when the sample tray assembly 112 has been loaded into the opening 124 and sample 
reading is to begin. Further details of the stage assembly 110 are described in the above-referenced U.S. Patent No. 
6,043.880 

[0028] To read the samples contained in the microwells 118 of a sample tray assembly 112 that has been loaded 
55 into the stage assembly 110, the stage assembly 110 is conveyed past a light sensing bar 130 as shown in Fig. 4. The 
light sensor bar 130 includes a plurality of light emitting/detecting ports 132. The light emitting/detecting ports 132 are 
controlled to emit light towards a column of 8 microwells 118 when the stage assembly 110 positions those microwells 
118 over the light emitting/detecting ports, and to detect fluorescent light being emitted from the samples contained in 
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those microwelis 118. In this example, the light sensor bar 130 includes 8 light emitting/detecting ports 132 which are 
arranged to substantially align with the 8 microwelis 118 in a column of the microwell array 116 when that column of 

microwelis 118 is positioned over the light emitting/detecting ports 132. 

[0029] The light emitting/detecting ports 132 are coupled by respective fiber optic cables 134 to respective light 
5 emitting devices 136, such as LEDs or the like. The light emitting/detecting ports 132 are further coupled by respective 
fiber optic cables 138 to an optical detector 140, such as a photomuttlplier tube or the like. Further details of the light 
sensor bar 130 and related components, as well as the manner in which the stage assembly 110 is conveyed past the 
light sensor bar 130 for reading the samples contained in the microwelis 118. are described in the above-referenced 
U.S. Patent No. 6,043.880. 

10 [0030] In general, one reading for each microwell is taken at a particular interval in time, and additional readings of 
each microwell are taken at respective intervals in time for a predetermined duration of time. In this example, one 
microwell reading is obtained for each microwell 118 at approximately one minute intervals for a period of one hour. 
One reading of each of the normalization wells 127, as well as one "dark" reading for each of the light emitting/detecting 
ports 132, are taken at each one minute intervals. Accordingly, 60 microwell readings of each microwell 118, as well 

15 as 60 readings of each normalization well 127 and 60 dark readings, are obtained during the one hour period. As 
discussed above, a reading is a measurement of the intensity of the fluorescent emission being generated by a micro- 
well sample in response to excitation light emitted onto the sample. These intensity values are stored in magnitudes 
of relative fluorescent units (RFU). A reading of a sample having a high magnitudeof fluorescent emissions will provide 
an RFU value much higher then that provided by a reading taken of a sample having low fluorescent emissions. 

20 [0031] Once the total number of readings (e.g. 60 readings) for each sample well have been taken, the readings for 
each sample must be interpreted by the well reading apparatus 100 so the well reading apparatus 100 can indicate 
whether the sample has tested positive or negative for a particular disease (e.g. CT or GC). The microprocessing unit 
of the well reading apparatus 1 00 is controlled by software to perform the following operations on the data representing 
the sample well readings. The operations being described are applied in essentially the same manner to the readings 

25 taken for each sample microwell 118. Accordingly, for illustrative purposes, the operations will described with regard 
to readings taken for one sample microwell 118, which will be referred to as the first sample microwell 118. 
[0032] As discussed above, during each one-minute interval in which all of the microwelis 118 in the sample tray 
assembly 112 are read, the light sensor bar 130 reads the normalization wells one time. Hence after 60 readings of 
each microwell sample have been taken, each normalization well 127 has been read 60 times by its respective light 

30 emitting/detecting port 132 of the light sensor bar 130. which results in eight sets of 60 normalization well readings. 
For illustrative purposes, the normalization readings of the normalization well 127 that has been read by the light 
emitting/detecting port 132 which has also read the first sample microwell 118 now being discussed, are represented 
as n^ through Hqq. 

[0033] Additionally, as discussed above, during each one-minute interval, the optical detector 140 is controlled to 
35 obtain a "dark" reading in which a reading is taken without any of the light emitting devices 136 being activated. This 
allows the optical detector 140 to detect any ambient light that may be present in the system. The dark readings are 
taken for each light emitting/detecting port 132. Accordingly, after 60 readings of every microwell 118 have been ob- 
tained, eight sets of 60 dark readings (i.e.. one set of 60 dark readings for each of the eight light emitting/detecting 
portions 132) have been obtained. For illustrative purposes, the dark readings obtained by the light emitting/detecting 
40 port 132 which read the first sample microwell 118 now being discussed are represented as through dgg- 

[0034] Fig. 5 is a graph showing the relationship of the 60 readings for one well which have been obtained during 
the one hour reading period. For illustrative purposes, these readings are represented as r^ through rgo- These readings 
are plotted on the graph of Fig. 5 with their RFU value being represented on the vertical axis with respect to the time 
in minutes at which the readings were taken during the reading period. 
45 [0035] As can be appreciated from the graph, the RFU values for the readings taken later in the reading period are 
greater than the RFU values of the readings taken at the beginning of the reading. For illustrative purposes, this example 
shows the trend in readings for a well that contains the particular disease (e.g.. CT or GC) for which the well is being 
tested. 

[0036] As can also be appreciated from Fig. 5, the graph of the "raw data" readings include a noise spike and a step 
50 as shown. The process that will now be described eliminates any noise spikes, steps or other apparent abnormalities 
In the graph which are the result of erroneous readings being taken of the sample well. 

[0037] The flowchart shown in Fig. 6 represents the overall process for interpreting the graph of raw data readings 
ri through rgQ shown in Fig. 5 to provide a well sample result which is used to determine whether the well sample 
includes the particular target disease for which it is being tested. These processes are performed by the controller (not 
55 shown) of the well reading apparatus 100 as controlled by software, which can be stored in a memory (not shown) 
resident in the well reading apparatus 100, or on a disk inserted into disk drive 106. 

[0038] As shown in Fig. 6, the software initially controls the controller to perform a dark correction on the normalizer 
data readings n^ through n^Q and on the well readings r-t through r^Q. The details of this step are shown in the flowchart 
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of Fig. 7 and in Step 1 of the pseudo-code set forth in the attached Appendix. 

[0039] In particular, in Step 1010, the dark reading values d^ through 6qq are subtracted from the corresponding 
normalizer reading values n^ through n^Q, respectively, to provide corrected normalizer readings cn^ through cngQ, 
respectively. That is, dark reading d-i is subtracted from normalizer reading n^ to provide corrected normalizer reading 

5 cn^, dark reading 62 is subtracted from normalizer reading ^2 to provide corrected normalizer reading cn2, and so on. 
[0040] The processing then proceeds to Step 1020 in which the dark readings d^ through 6qq are subtracted from 
their corresponding well readings r-, through reo. respectively to provide corrected well readings c^ through Cqq, respec- 
tively. That is. dark well reading d^ is subtracted from well reading r^ to provide corrected well reading c-,, dark reading 
62 is subtracted from well reading r2 to provide corrected well reading cr2, respectively, and so on. 

10 [0041] After all of the corrected normalizer readings and corrected well readings have been obtained, the processing 
continues to the filtering operations Step 1100 of the flowchart shown in Fig. 6, in which noise is filtered from the 
corrected normalizer readings cn^ through c^qq, which were obtained during Step 1010 described above. Details of 
Step 1100 are shown in Fig. 8. and in Step 2 in the attached pseudo-code. Specifically, in this example, a 5 point 
running median is applied to the corrected normalizer readings cr^ through ctqq. 

15 [0042] As shown in Step 1110, the first two smoothed nomrializer values xn^ and xn2 are set equal to the first two 
corrected normalizer values cn^ and cn2, respectively, while the last two smoothed normalizer values xn^Q and xngo 
are said equal to the last two corrected normalizer values cn^g and cngQ. Then, In Step 1120, smoothed normalizer 
values xng through xn53 are obtained as an average or midpoint value of their corresponding corrected normalizer 
values cn3 through cn^, respectively, and surrounding corrected normalizer values. For example, smoothed normalizer 

20 value xn3 is set equal to an average or midpoint value of corrected normalizer values cn^, cng, cn3, cn4, and eng. 
Similarly, the smoothed normalizer value xn4 is set equal to an average or midpoint value of corrected normalizer 
values eng. cn3, cn4, cng, and cng, and smoothed normalizer values xng through xn^Q are calculated in a similar manner. 
[0043] Once all smoothed normalizer values xn^ through xngo have been obtained, the processing continues to the 
dynamic normalization step 1200 shown in the flowchart of Fig. 6. The details of the dynamic normalization process 

25 is shown in the flowchart of Fig. 9, as well as in Step 3 of the attached pseudo-code. Specifically in this example, the 
smoothed normalizer values xn-, through xngQ, as well as the corrected well reading values cr^ through ergo are used 
to calculate dynamic normalization values in nr^ through nrgQ. 

[0044] In Step 1210, a scalar value is set which is employed in the calculations. In this example, the scalar value is 
3000, but can be any suitable value. The processing then proceeds to Step 1220. where the scalar value, corrected 

30 well reading values, and smoothed normalize values are used to calculate dynamic nomrialization values. In particular, 
to calculate the dynamic normalization values, the corresponding corrected well value is multiplied by the scalar value, 
and then that total is divided by the corresponding smoothed normalizer value. For instance, to obtain dynamic nor- 
malization value nr^, corrected well reading value cr^ is multiplied by 3000 (the scalar value), and then that total is 
divided by the value of smoothed normalizer xn^ . Similarly, dynamic normalization value nr2 is calculated by multiplying 

35 corrected well reading value cr2 by 3000, and then dividing that total by smoothed normalizer value xn2. This process 
continues until all 60 dynamic normalization values nr^ through nrgo have been obtained. 

[0045] The processing then continues to perform the input noise filtering operation on the well data as shown in Step 
1300 of the flowchart in Fig. 6. The details of this operation are shown In the flowchart of Fig. 10. as well as in Step 4 
of the attached pseudo-code. In Step 1300, a three point running median is applied to the dynamic normalization values 

40 nr^ through nrgo to obtain smoothed normalized values x^ through Xqq. To perform this operation, as shown in Step 
1310. the first smoothed normalized value x^ is set equal to the first dynamic normalization value in r^, and the last 
smoothed normalized value x^q is set equal to the last dynamic normalization value in r^Q. The processing then proceeds 
to Step 1320, where the smoothed normalized values X2 through X59 are obtained. These values are obtained by 
applying a three point running median to the dynamic normalization values, such that the smoothed normalized values 

45 X2 through X59 are obtained by calculating an average or midpoint of their corresponding dynamic normalization values 
nr2 through nrgg, respectively, and the surrounding normalization values. 

[0046] That is, in this example, smoothed normalized value X2 is obtained by taking the average or midpoint of 
dynamic normalization values nr^, nr2, and nr3. Similarly, smoothed normalized value X3 is obtained by taking the 
average or midpoint of dynamic normalized value nr2. nr^. and nr4. Smoothed normalize values X4 through X59 are 

50 obtained in a similar manner. 

[0047] Once the smoothed normalized values x-, through Xqq have been obtained, a three point running median is 
applied to those values to obtain smoothed normalized values z^ through z^q. That is, in Step 1330, smoothed nor- 
malized value z^ is set equal to smoothed normalized value x^, and smoothed normalized value z^q is set equal to 
smoothed normalized value Xqq, Then, in Step 1340, smoothed normalized values Z2 through Z59 are obtained by 

55 calculating an average or midpoint of their corresponding smoothed normalized value X2 through X59 and the surround- 
ing smoothed normalized values. That is, smoothed normalized value Z2 is obtained by calculating the average or 
midpoint of smoothed normalized values x-j, X2, and X3. Similarly, smoothed normalized value Z3 is obtained by calcu- 
lating the average or midpoint of smoothed normalized values X2, X3, and X4. Smoothed normalized value Z4 through 




7 



EP1 158 449 A2 

are then obtained in a similar manner. 
[0048] After Steps 1000 through 1300 of the flowchart in Fig. 6 have been performed as described above, the well 
readings have therefore been smoothened and normalized, and are represented by the smoothed normalized values 
Zi through Z50. Accordingly, as shown in the graph of Fig. 11, when the smoothed normalized values through Zqq 
5 are plotted with respect to a corresponding time periods in which their corresponding well readings have been obtained, 
the noise spikes in the graph have been eliminated. 

[0049] However, these smoothing and normalizing operation did not remove the step which is present in the graph 
as shown in Fig. 11. This increase in the reading values, which resulted in the step appearing in the graph, was likely 
caused by the presence of a bubble in the well which formed after the 30^^ well reading was obtained (i.e., after an 
10 elapsed time of 30 minutes), but before the 31 well reading was obtained. Accordingly, the magnitude of well reading 
values through r^g and hence, and magnitude of smoothed and normalized values z^^ through Zgo have been 
increased due to the presence of this bubble. Therefore, it is necessary to reduce the smoothed normalized values 
Z3^ through z^q by a value proportionate to the size of the step. 

[0050] The step removal operation is performed in Step 1400 as shown in the flowchart In Fig. 6. Details of the step 

15 removal operation are set forth in the flowchart in Fig. 12, and in Step 5 of the attached pseudo-code. 

[0051] It has been determined that graphs of these types generally will have only one or possibly two steps, and will 
almost never have more than five steps. Accordingly, all of the steps in the graph will have been located and removed 
after performing the step locating process five times. Accordingly, in Step 1 405 in the flowchart of Fig. 1 2, a count value 
is set to allow the process to repeat a maximum of five times. The process then proceeds to step 141 0, where difference 

20 values dr^ through dr^g are calculated which represent the differences adjacent smoothed normalized value z^ through 
ZgQ. That is, the first difference value dr^ is calculated as the value of smoothed normalized value Z2 minus smoothed 
normalized value z^. The second difference value dr2 is calculated as the value of smoothed normalized value Z3 minus 
smoothed normalized value Z2. This process is repeated until 59 difference values dr^ through drgg have been obtained. 
[0052] The processing then continues to Step 1415, in which the difference values dr^ through drgg are added together 

25 to provide an average total, which Is then divided by 59 to provide a difference average 'dr. The processing then 
continues to Step 1420, where a variance value var(dr) is calculated. This variance value is calculated by subtracting 
the difference value 'drfrom each difference value dr^ through drgg, squaring each subtraction value, and then summing 
the totals of the squared values. For example, the difference value 'dr is subtracted from the first difference value dr-| 
to provide a total, which is then squared. The difference value *dr is then subtracted from the second difference value 

30 dr2, and that total is squared. This process continues for all remaining difference values dr3 through drsg. These 59 
"squared" totals are then added and divided by 58 to obtain the variance value (dr). 

[0053] The process then continues to Step 1425 where a sum value "s" is calculated. This sum value is calculated 
by subtracting the difference value 'dr from each of the difference values dr^ through drsg, taking each result to the 
fourth power to obtain a set of 59 quadrupled results, and then adding all of the 59 quadrupled results. That is, the 

35 difference value 'dr is subtracted from the first difference value dr^ to provide a result. That result is then taken to the 
fourth power to provide a first quadrupled result. The difference value 'dr is subtracted from second difference value 
dr2, and the result of the subtraction is taken to the fourth power to provide a second quadrupled result. This process 
is repeated for the remaining difference values dr3 through drsg until all 59 quadrupled results have been calculated. 
The 59 quadrupled results are then added to provide the sum value "s". 

40 [0054] In Step 1430, the processing determines whether the variance value var(dr) is equal to zero. If the value of 
var(dr) is equal to zero, the processing proceeds to Step 1433, where the count value is incremented by one, and 
Steps 1410 through 1425 are repeated as discussed above. However, if the value of var(dr) is not equal to zero, then 
the processing proceeds to Step 1435. In Step 1435, as critical value CRIT_VAL is set equal to 4.9. The processing 
then proceeds to Step 1440, where it is determined whether the value of sum "s" divided by the total of var(d) squared 

45 multiplied by 59 is greater than the value of CRIT_VAL. If the calculated value is not greater than CR1T_VAL. then the 
step location and repair processing is completed, and the processing continues to the periodic noise filter processing 
in Step 1500 shown in the flowchart of Fig. 6. 

[0055] However, If the total is greater than the value of CRIT_VAL, then the processing proceeds to Step 1 445 where 
processing will be performed to determine the location of the step. This is accomplished by subtracting the difference 

50 value 'dr from each of the 1 through 59 difference values dr^ through 6r^Q, taking the absolute value of each of those 
subtraction results, and determining which of the absolute values is the greatest. For example, the processing first 
subtracts the values 'dr from the first difference value dr^, and takes the absolute value of that result. This absolute 
value is compared with a variable "max", which has Initially been set to zero. If the absolute value is greater, the variable 
max is set to that absolute value, and the variable maxpt_dr is set equal to the number of the difference value, which 

55 in this case is 1 . 

[0056] The processing then subtracts difference value 'dr from the second difference value dr2, and takes the absolute 
value of that subtraction value. The processing determines whether that absolute value is greater than the new "max" 
value. If the absolute value is greater, then "max" is set to that absolute value, and maxpt_dr is set equal to 2. This 
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process is repeated for all remaining difference values 6t^ through drgg. After the process Is completed. maxpt_dr will 
be set equal to the number of the smoothed normalized value at which the step has occurred. As discussed above, in 
this example, it is presumed that the step occurred at value 230. Accordingly, maxpt_dr is set to 30. 
[0057] The process then continues to Step 1450 during which the median difference value of the difference values 

5 dr^ through drgg is determined. Then, in Step 1455, the smoothed normalized values occurring after the step are de- 
creased by the difference value calculated for the smoothed normalized value at which the step occurred, and then 
Increased by the median difference value calculated in Step 1450. For example, the smoothed normalized values z^^ 
through Zqq are each decreased by the magnitude of difference drgQ (the step occurred after the 30^*^ reading), and 
then the smoothed normalized values 23^ through z^q are then each increased by the median difference value calculated 

10 in Step 1450. As shown in Fig. 13, this process has the affect of shifting the entire portion of the curve representing 
the RFU values of 23^ through z^q downward, thus eliminating the step. 

[0058] The processing then proceeds to Step 1460 where It is determined whether the entire process has been 
repeated five times. If the value of count does not equal five, the value of count is increased by one in Step in 1465, 
and the processing returns to Step 1410 and repeats as discussed above. However, if the value of count is equal to 
15 five, the processing proceeds to the periodic noise filter Step 1500 In the flowchart shown In Fig. 6. 

[0059] The periodic noise filtering operation 1 500 is performed to further filter out erroneous values which may exist 
in the graph shown in Fig. 13 in which the step has been repaired. Details of the periodic noise filtering operation are 
shown in the flowchart of Fig. 14 and in Step 6 of the attached pseudo-code. 

[0060] Specifically, a five-point moving average is applied to the read values z^ through 250 represented in the graph 
20 of Fig. 13 to provide filtered values f^ through fgo. In Step 1510, the first two filter values f-, and fg are set equal to the 
smoothed normalized values z^ and Z2, respectively, and the last two filtered values f^g and f^Q are set equal to the 
two smoothed normalized values 259 and Zgo, respectively. Then, in Step 1520, the filtered value fg through fgg are 
determined by taking the average of corresponding smoothed normalized values Z3 through Z58, respectively and 
surrounding smoothed normalized values. For example, filtered value is determined by taking the sum of smoothed 
25 normalized values z^, Z2, 23, Z4, and 25, and dividing that sum by 5. Filtered value ^4 is determined by taking the sum 
of smoothed normali2ed value 23, 23, 24, 25, and 2g, and dividing that sum by 5. This process Is repeated until all 
remaining filter f3 through fgg have been obtained. 

[0061] The processing then continues to Step 1600 shown in Fig. 6, in which the processing determines whether 
the filtered values f^ through fgQ. which were derived from the above-described steps from the raw data well read values 
30 r-i through t^q, respectively, were actually taken of a well, or, In other words, whether a well was actually present at that 
location in the microwell array 116 of the sample tray assembly 112. Detailsof the well present determination processing 
are shown in the flowchart of Fig. 15 and in Step 7 of the pseudo-code. 

[0062] Specifically, in Step 1610, a well reading average wpgyg is determined by adding the filter values f^Q. f2o. ^30' 
f40 and fgQ, and dividing those values by 5. This well present average wpgyg is compared to a well threshold value 

35 WP_THRES, which in this example is set to 125.0. If. in Step 1620, the processing determines that the well present 
average wpgyg is greater than 2ero and less than the threshold value WP_THRES, then the processing determines 
that no well is present and that the data obtained is entirely erroneous. The processing then proceeds to Step 1900 in 
the flowchart shown in Fig. 6 where processing for that well is concluded. However, if the processing determines in 
Step 1620 that a well is present, the processing continues to Step 1700 In the flowchart shown in Fig. 6. 

40 [0063] In Step 1 700, the processing establishes a base line background correction, in which an average value based 
on, for example, the first five background values f^ through fg, is calculated. Other ranges of values, such as f^o through 
fi5, may be used, depending on the assay. This average value is then subtracted from all of the filtered values f^ through 
fgo- Additionally, the values used to calculate the average value can each be set to zero after being used to calculate 
the average value, although this is not required. Further details of this background correction operation Is shown In the 

45 flowchart of Fig. 16 and in Step 8 of the pseudo-code. 

[0064] That is. in Step 1710. the filtered values fi through fg are added to produce a sum value. In this example, the 
sum value is then divided by 5 to provide an initial adjustment value lA. 

[0065] The processing then proceeds to Step 1720, in which the initial adjustment value lA is subtracted from each 
of the filter values f^ through fgo- If, in performing the subtraction, the filter value becomes less than zero, that filtered 
50 value is set to zero. As shown in the graph of Fig. 17, this processing shifts the portion of the graph between filter 
values fi and Iqq down toward the horizontal axis. The processing then proceeds to Step 1800 in the flowchart of Fig. 
6 to compare the corrected graph to a threshold value as will now be described. 

[0066] It Is noted that all of the processing steps performed above remove extraneous and additive values to provide 
a graph most suitable for comparison to the threshold value. However, it is only necessary to perform the background 
55 correction to produce a graph which can be compared to the threshold value to provide acceptable results. 

[0067] The threshold comparison step is described in more detail In the flowchart shown in Fig. 18. In step 1810; the 
corrected graph is compared to a threshold value which has been determined based on past results. In other words, 
it has been determined from past data readings taken to detect the presence of a particular disease (e.g.. CT or GC) 
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that the corrected reading values will exceed a particular value at some time during the 60 readings if that disease is 
indeed present in the sample (i.e., the sample is positive). Accordingly, the threshold value is chosen to be a value 
which will provide the most accurate indication as to the presence or absence of a disease of interest in the sample. 
This can be accomplished by choosing a threshold value which simultaneously maximizes sensitivity and specificity. 

5 In this example, the threshold value is chosen to be "100". 

[0068] If the comparison determines that the first corrected value to exceed the threshold value is a value other than 
the first several values used to calculate the background correction value, the processing proceeds to steps 1820 and 
1830 as shown in Fig. 18. The processing then interpolates to determine where the graph crosses the threshold. That 
is, as shown in Fig. 19, the value at which the graph crosses the threshold is between two actual reading values. 

10 Accordingly, the processing interpolates between the two actual values to calculate a reported value which is a value 
corresponding to the position on the x-axis (the time axis). For instance, if the graph crosses the threshold between 
readings 1 6 and 1 7, the interpolation will arrive at a reported value having a magnitude between 1 6 and 1 7 (e.g., 1 6.40). 
The processing then proceeds to step 1840, where the controller controls the well reading apparatus 100 to report the 
reported value and provide an indication that the sample in the corresponding well has tested positive for the targeted 

15 disease. This indication can be in the form of a display on the display screen 108, in the form data stored to a disk in 
the disk drive 106, and/or in the form of a print-out by a printer resident in or attached to the well reading apparatus 100. 
[0069] However, if the processing determines in step 1810 that the first corrected value to exceed the threshold value 
is a value among the first several values used to calculate the background correction value, the processing proceeds 
to steps 1850 and 1860 as shown in Fig. 18. This condition is shown in the graph of Fig. 20, and occurs when the 

20 sample is a very high positive sample containing a substantial amount of the pathogen of interest or target nucleic 
acid. The processing then interpolates between the last value among the values used for background correction and 
the next value outside the background correction value to calculate a reported value. For example, if the last value 
among the values used for background correction is 5, the interpolation will calculate a value between 5 and 6 (e.g., 
5.70). The processing then proceeds to step 1840, where the controller controls the well reading apparatus 100 to 

25 provide an indication that the sample in the corresponding well has tested highly positive for the targeted pathogen or 
nucleic acid, in a manner similar to that described above. This indication can be in the form of a display on the display 
screen 108, in the form data stored to a disk in the disk drive 106, and/or in the form of a print-out by a printer resident 
in or attached to the well reading apparatus 100. 

[0070] If the processing determines in step 1810 that none of the values in the graph exceed the threshold value, 
30 the processing proceeds to steps 1870 and 1880 as shown in Fig. 18. This condition is shown in the graph of Fig. 21, 
and occurs when the sample is a negative sample containing none or an insufficient amount of the pathogen of interest 
or target nucleic acid. The processing then designates the last value in the graph (i.e., 60) as the reported value. The 
processing then proceeds to step 1840, where the controller controls the well reading apparatus 100 to provide an 
indication that the sample in the corresponding well has tested negative for the targeted disease, in a manner similar 
35 to that described above. This indication can be in the form of a display on the display screen 108, in the form data 
stored to a disk in the disk drive 106, and/or in the form of a print-out by a printer resident in or attached to the well 
reading apparatus 100. 

[0071] As discussed above, the manner in which the samples from patient number 1 collected in the other sample 
microwells are read and analyzed is essentially identical to that described above for the sample in the first sample 
40 microwell. Specifically, the 60 readings taken of the sample in each of the respective sample microwells are processed 
according to Steps 1000 through 1800 in Figs. 6 as described above. 

[0072] The above processing can then performed for all of the patient samples in essentially the same manner. As 
discussed above, if each patient sample is being tested for 2 diseases, the microwell array 116 and the sample tray 
assembly 112 can accommodate samples from a maximum number of 30 patients. However, in some instances, the 

45 amplification control AC can be included as an internal control present in all of the wells. This amplification control can 
be illuminated by a different frequency light emitted by the tight emitting/detecting ports 132 of the light sensor bar 130, 
or by a duplicate light sensor bar (not shown) as described in more detail in, for example, the above-referenced U.S. 
Patent No 6.043,880. In this event, a set of two microwells for each patient is required to test for two different diseases 
(e.g., CT and GC), and only one microwell per patient is required to test for one disease (e.g., CT or AC). 

50 [0073] It is also noted that before any results are reported to patients, the values obtained from reading the CT, GC 
and AC positive and negative control samples are processed in the manner described above with regard to Steps 1000 
through 1800. and the resulting values are analyzed to assure that the known positive and negative control samples 
have indeed been read as positive and negative samples, respectively. If the readings of any of these control samples 
are incorrect (I.e., a negative sample has been identified as a positive sample or vice-versa), all of the sample readings 

55 taken for the entire microwell tray are called into question. All of the sample data must be discarded, and new samples 
can be gathered in a new microwell array, and then read and evaluated in the manner described above. 
[0074] Although only a few exemplary embodiments of this invention have been described in detail above, those 
skilled in the art will readily appreciate that many modifications are possible in the exemplary embodiments without 
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materially departing from the novel teachings and advantages of this invention. Accordingly, all such modifications are 
intended to be included within the scope of this invention as defined in the following claims. 
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APPENDIX 



Pseudo-code 

10 



RFQ - (r^ . . rco) for one well 
IStBFO ■ (Oi . . n«g) for nozmalizor column 
15 DRFO - (dx . . d«a) dark for entire plate 

1. Dark Correct Readings 

for i - 1 to €0 

// Dark correct normal izer readings 
cni - ni - di 

20 

// Dark correct well readings 
cri « ri - di 

next i 

2. Smooth Normal! zers 

25 ^ply 5 point running median to dark corrected normalizers 

for i • 1 to 56 

// Sort en* , cn^^j , cn^^, , c%^, , tn^^ and Set xa^^^ « c^^nia^ 

next i 

30 3. Dynamic Normalization 

scalar - 3000 

for 1 « 1 to 60 

nr^ « (cr^ * scalar) / xn^ 

next i 

35 

4. Smooth normalised data 

Apply 3 point running median to normalized data 
^ - nr^,. x«o - nr,o 
for i • 1 to 58 

// Sort nr^ , nr^^^ , nr^^, and Set Xj^^ • nr, 

40 next 1 

Apply 3 point running median to normalized data again 
for i « 1 to 58 



45 



50 



55 



next 



^ // Sort Xi , xi^i , Xj^a and Set z^^^ « 



5. Step location and repair 
for count * 1 to S 
a) Sat up diffarences between adjacent points, 
for i • 1 to S9 

next i^'^ - - 

iim''!^r^*'=^ differences between points, 

for i - 1 to 59 

sum - sum dr. 
next i ' 
*dr - sum / 59 
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c) Calculate variance 
su;n • 0 

for i - 1 to 59 

sum • a\im +• (dr^ - *dr)* 

next i 

var(dr) • aum / 58 Note: 58 

aum « 0 

for i - i to 59 

aum ■ Bum ♦ (dr^ - *dr)* 

next i 
a ■ aum 

if var(dr) - 0 then next coxint 
CRZT VAL - 4*9 

if a""/ ( 59 * varCdr)* ) > CRVTjmu then 
// Step Occurred, Iiocate where 
max ■ 0 

for i - 1 to 59 - . . 

7/ Note: Location method inaurea finding the I"' maximum 
// difference 
if jdr^ - *dr| > max then 
max « |dri - "dr| 
* ^ raaxpt_dr « i 

•end if 

next i 

Replace max dr point by median difference: 

Sort dr^ , dr, , . . dr„ and Set dr^^ 

// From atep point onward, adjuat curve down by dr,„^_a. and then 
// back up by dr,^ 

for i - (maxpt^dr-fl) to 60 Note: maxpt_dr+l 

■ - dr^^ ^ dr. 

next i 



40 



45 



50 




maxpt_dr' ' maxpt_di+l 



^lee 

Break out of count loop 

end if 
next count 



55 



13 



15 



20 



25 



EP 1 158 449 A2 



6. Filter 

;^ply S point moving average to repaired data 
* for i « 1 to 56 

next i 

7. Apply Well Present Threshold (Re£. 9) 

10 wp^ - (£ia 4. £„ ♦ ♦ + ffo) / 5 

HPJTKRES - 125.0 

if wp^ > 0 and vp^ < HPJXHRBS then 
// No well is present 
Exit calculation 

end if 



8. Establish baseline 

Average RFC values from to K, « where > and are parameters 
currently set to be l: and 5^ 

Initial adjustment 

// Ayg points^ where Kj ■ 1 and K, ■ 5" 
sum >• 0 
for i M to 

pURi m sun * f £ 

next i 
XA m sum / (K, - 1) 

// 
// 

30 for i « <Ks > 1) to 60 // ,/l to 60 

f i - fi • lA 
if f^ < 0 then f A - 0 

next i 
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Claims 



40 



1 . A computerized method for controlling a system to analyze numerical data pertaining to a sample assay comprising 
at least one biological or chemical sample, and numerical data including a set of data pertaining to each respective 
sample, each said set of data including a plurality of data values, each representing a condition of said respective 
sample read at a respective point in time, said method comprising the steps of: 

for each sard set of data, performing the steps of: 

assigning a respective numericai value to each of said data values; 
calculating a correction value based on at least one of said numerical values; 

adjusting each of said numerical values based on said correction value to provide respective corrected nu- 
50 merical values; 

determining whether any of said corrected numerical values exceeds a threshold value; 

and 

controlling said system to Indicate v\/hether said sample has a predetermined characteristic based on a result 
55 of said determining step. 

2. A method as claimed in claim 1. wherein: 
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said assigning step includes the step of arranging said data values in a sequence representative of said re- 
spective points in time; and 

said determining steps sequentially compares said corrected numerical values to said threshold value until 
one of said corrected numerical values exceeds said threshold value. 

5 

3. A computerized method for controlling a system to analyze numerical data pertaining to a sample assay comprising 
at least one biological or chemical sample, said numerical data including a set of data pertaining to each respective 
sample, each said set of data including a plurality of data values, each representing a condition of said respective 
sample read at a respective point in time, said method comprising the steps of: 

10 for each said set of data, performing the steps of: 

representing each of said plurality of data values as points on a graph having a vertical axis representing the 
magnitudes of said values and a horizontal axis 

representing a period of time during which readings of said sample were taken to obtain said plurality of data 
15 values; 

calculating a correction value based on at least one said magnitudes of said data values; 

adjusting each of said data values based on said correction value to provide a corrected graph of respective 

corrected value; 

determining whether any of said corrected values of said corrected graph exceeds a threshold value in a 
20 direction of said vertical axis; and 

controlling said system to indicate whether said sample has a predetermined characteristic based on a result 
of said determining step. 

4. A method as claimed in claim 3, wherein: 

25 said calculating step calculates said correction value based on a plurality of said magnitudes of said data 

values at a beginning of said graph. 

5. A computer-readable medium of Instructions for controlling a system to analyze numerical data pertaining to a 
sample assay comprising at least one biological or chemical sample, said numerical data including a set of data 

30 pertaining to each respective sample, each said set of data including a plurality of data values, each representing 

a condition of said respective sample read at a respective point In time, said medium of instructions comprising: 

a first set of instructions, adapted to control said system to assign a respective numerical value to each of said 

data values; 

35 a second set of instructions, adapted to control said system to calculate a correction value based on at least 

on of said numerical values; 

a third set of instructions, adapted to control said system to adjust each of said numerical values based on 
said correction value to provide respective corrected numerical values; 

a fourth set of instructions, adapted to control said system to determine whether any of said corrected numerical 
40 values exceeds a threshold value; and 

a fifth set of instructions, adapted to control said system to indicate whether said sample has a predetermined 
characteristic based on a result of said determining step. 

6. A computer-readable medium of instructions as claimed in claim 5, wherein: 

45 

said first set of instructions is further adapted to control said system to arrange said data values in a sequence 
representative of said respective points in time; and 

said fourth set of instructions is further adapted to control said system to sequentially compare said corrected 
numerical values to said threshold value until one of said corrected numerical values exceeds said threshold 
50 value. 

7. A system for analyzing numerical data pertaining to a sample assay comprising at least one biological or chemical 
sample, said numerical data including a set of data pertaining to each respective sample, each said set of data 
including a plurality of data values, each representing a condition of said respective sample read at a respective 

55 point in time, said system comprising: 

means for assigning a respective numerical value to each of said data values; 

means for calculating a correction value based on at least one of said numerical values; 
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means for adjusting each of said numerical values based on said correction value to provide respective cor- 
rected numerical values; 

means for determining whether any of said corrected numerical values exceeds a threshold value; and 
means for controlling said system to indicate whether said sample has a predetermined characteristic based 
5 on a result of said determining step. 



8. A system as claimed in claim 7. wherein: 



said assigning means includes means for arranging said data values in a sequence representative of said 
10 respective points in time; and 

said determining means sequentially compares said corrected numerical values to said threshold value until 
one of said corrected numerical values exceeds said threshold value. 



9. A system for analyzing numerical data pertaining to a sample assay comprising at least one biological or chemical 
sample, said numerical data including a set of data pertaining to each respective sample, each said set of data 
including a plurality of data values, each representing a condition of said respective sample read at a respective 
point in time, said system comprising: 



means for representing each of said plurality of data values as points on a graph having a vertical axis repre- 
senting the magnitudes of said values and a horizontal axis representing a period of time during which readings 
of said sample were taken to obtain said plurality of data values; 

means for calculating a correction value based on at least one said magnitudes of said data values; 
means for adjusting each of said data values based on said correction value to provide a conrected graph of 
respective corrected values; 

means for determining whether any of said corrected values of said corrected graph exceeds a threshold value 
in a direction of said vertical axis; and 

means for controlling said system to indicate whether said sample has a predetermined characteristic based 
on a result of said determining step. 



10. A system as claimed in claim 9, wherein: 

said calculating means calculates said correction value based on a plurality of said magnitudes of said data 
values at a beginning of said graph. 
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